The physical determinants of the thickness of lamellar polymer crystals 
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Based upon kinetic Monte Carlo simulations of crystallization in a simple polymer model we 
present a new picture of the mechanism by which the thickness of lamellar polymer crystals is 
constrained to a value close to the minimum thermodynamically stable thickness. This description 
contrasts with those given by the two dominant theoretical approaches. 
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On crystallization from solution and the melt many 
polymers form lamellae where the polymer chain tra- 
verses the thin dimension of the crystal many times fold- 
ing back on itself at each surface jj]. (The crystal ge- 
ometry is shown by the example configuration in Figure 
|l|). Although the lamellar crystals were first observed 
over forty years ago their physical origin is still contro- 
versial. It is agreed that the kinetics of crystallization 
are crucial since extended-chain crystals are thermody- 
namically more stable than lamellae. However, the ex- 
planations for the dependence of the lamellar thickness 
on temperature offered by the two dominant theoretical 
approaches appear irreconcilable || . The lamellar thick- 
ness is always slightly greater than Z m i n , the minimum 
thickness for which the crystal is thermodynamically sta- 
ble; Z m j n is approximately inversely proportional to the 
supercooling. 

The first theory, which was formulated by Lauritzen 
and Hoffman soon after the initial discovery of the chain- 
folded crystals invokes surface nucleation of a new 
layer on the thin side faces of the lamellae as the key 
process. It assumes that there is an ensemble of crys- 
tals of different thickness each of which grows with con- 
stant thickness; the thickness for which the crystal has 
the maximum growth rate is then sought. The growth 
rates are derived by assuming that a new crystalline layer 
grows by the deposition of a succession of stems (straight 
portions of the polymer that traverse the crystal once) 
along the growth face. The two main factors that de- 
termine the growth rate are the thermodynamic driving 
force and the free energy barrier to deposition of the first 
stem in a layer. The former only favours crystallization 
when the thickness is greater than l m { n ; the latter in- 
creases with the thickness of the crystal because of the 
free energetic cost of creating the two new lateral surfaces 
on either side of the stem and makes crystallization of 
thick crystals prohibitively slow. Therefore, the growth 
rate passes through a maximum at an intermediate value 
of the thickness which is slightly greater than ^ m i n - 

The second approach which has been termed the en- 
tropic barrier model, is based upon the interpretation of 
kinetic Monte Carlo simulations BM and rate-theory cal- 



culations [|6| of a simplified model of the polymer crystal 
growth. As with the surface nucleation approach, the ob- 
served thickness is suggested to result from the competi- 
tion between a driving force and a free energy barrier con- 
tribution to the growth rate. However, a different cause 
for the free energy barrier is postulated. As the polymer 
surface in the model can be rough, it is concluded the 
details of surface nucleation of new layers are not impor- 
tant. Instead, the outer layer of the crystal is found to 
be thinner than in the bulk; this rounded crystal profile 
prevents further crystallization. Therefore, growth of a 
new layer can only begin once a fluctuation occurs to an 
entropically unlikely configuration in which the crystal 
profile is 'squared-off '. As this fluctuation becomes more 
unlikely with increasing crystal thickness, the entropic 
barrier to crystallization increases with thickness. 

Although both approaches are able to describe cor- 
rectly some of the basic phenomenology of polymer crys- 
tallization, both have questionable aspects. For example, 
the surface nucleation assumption that crystals grow 




FIG. 1. Cut through a polymer crystal which was pro- 
duced by the growth of twenty successive layers on a sur- 
face with a uniform thickness of 50 units at a temperature, 
T = 2.0 efc -1 . Folds occur at the top and bottom surfaces, 
and the stems are represented by vertical cuboids. The cut 
is 16 stems wide. 
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with constant thickness is contradicted by experiments in 
which temperature changes cause the thickness of grow- 
ing lamellae to adjust to the new temperature 0. And in 
the entropic barrier approach it is not clear whether ap- 
proximations such as the implicit representation of the 
chain connectivity and chain folds by a set of simple 
growth rules compromise its conclusions. 

The approach we take combines some of the more at- 
tractive features of the two models. As in the surface 
nucleation model, we form new polymer layers by the 
sequential growth of stems, and growth is considered to 
proceed layer by layer. However, we remove any con- 
straints on the stem length. This additional complexity 
necessitates the use of simulations to characterize the re- 
sulting behaviour. Others have recognized the need for 
such a 'multiple pathway' approach, but, because of lim- 
ited computational resources at that time, only approxi- 
mate treatments were possible ||-[ll]] . 

We model the polymer by a self-avoiding walk on a 
simple cubic lattice. There is an attractive energy, -e, 
between non-bonded polymer units on adjacent lattice 
sites and between the polymer and the surface, and an 
energetic penalty, e g , for kinks ('gauche bonds') in the 
chain. e g defines the stiffness of the chain. Here, we use 
e g = 8e, however similar results are obtained at any pos- 
itive e g . We only model the crystalline portion of the 
polymer explicitly and the rest of the chain is assumed 
to behave like a two-dimensional ideal coil adsorbed onto 
the surface [ ^2[ . 

At each step in our kinetic Monte Carlo simulation 
p3[ changes of configuration can only occur at the ends 
of the crystalline portion of the polymer. The possible 
processes are for the crystalline layer to grow by one unit 
by the extension of one of the stems or by the creation of 
a fold, or for the crystalline part of the polymer to shrink 
by one unit. The free energy change for a step is given 
by A A = A_E xt ai + A^4 co ii, where Ai^tai is the change 
in energy of the crystalline configuration and AA co ii is 
the change in free energy of the coil. From this we derive 
a rate for a transition between the current state i and 
the state j: kij = min (1, exp (— AAij/kT)). We then 
randomly choose one of the states connected to i, with 
probabilities given by pij = kij/ Y]j, kij> and update the 
time by an increment At = — log(p) / ., kij' , where p is 
a random number in the range [0,1]. 

In this way we can grow a new crystalline layer on a 
surface which represents the growth face of the polymer 
crystal. By taking averages over many such layers we can 
statistically characterize the properties of the new layer. 

In Figure ^| we show how the thickness of the new layer 
depends both on temperature and the thickness of the 
underlying layer. We immediately see that the thickness 
of the new layer is not necessarily the same as that for 
the previous layer, showing that the constant thickness 
assumption of the surface nucleation approach does not 
hold for our model. For example, for growth on a surface 



that is 100 units thick, the new layer is always thinner. 
As the temperature approaches the melting temperature 
(T m = 4.06 efc -1 ) the thickness rises due to the increase 
of imin- At low temperature the thickness also increases 
because it becomes harder to scale the free energy barrier 
for forming a fold, and so on average the stem continues 
to grow for longer. 
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FIG. 2. The average stem length in the new crystalline 
layer as a function of temperature. The different curves are 
for different thicknesses of the underlying surface, as labelled. 

A different perspective is obtained by plotting the 
thickness of the new layer against the thickness of the 
underlying surface, as in Figure |^. It also allows us to 
consider multi-layer growth implicitly. At T = 2.75 efc -1 
growth on a 50-unit thick surface produces a new layer 
of thickness ~36. Growth of another layer on top of this 
layer is similar to growing a layer on a surface with a 
uniform thickness of 36; the latter gives a new layer of 
thickness ~28. We can see the effect of adding each suc- 
cessive layer by following the dotted line in Figure ||a; 
it also clearly shows that the mapping is a fixed-point 
attractor. The thickness of the crystal converges to the 
thickness at which the curve crosses y = x, i.e. to the 
point where the thickness of the new layer is the same 
as the previous. Therefore, at this temperature, there is 
only one thickness for which stable growth occurs (the 
inset of Figure |^a shows that this is not the thickness at 
which growth of a new layer is the most rapid) and there 
is no barrier to growth due to a rounded crystal profile. 

Although our results contrast with the assumptions of 
both surface nucleation and entropic barrier models that 
the observed thickness corresponds to the one with the 
maximum growth rate, in a number of the early polymer 
crystallization papers p|Js|JT4]| it was realized that stable 
growth could only occur at the thickness for which a new 
layer has the same thickness as the previous. However, 
since then this insight has for the most part been ne- 
glected, and, to the best of our knowledge, the iterative 
maps that underlie it have not been previously visualized. 
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FIG. 3. The average stem length in the new crystalline 
layer as a function of the thickness of the underlying surface 
for (a) T = 2.75 efc~\ (b) T = 3.375 e/T 1 . The dotted lines 
show how the thickness changes on addition of successive 
layers to a 50-unit thick surface. The inset of (a) shows the 
growth rate (polymer units per unit time) of the new layer 
as a function of thickness for T = 2.75 efc -1 . 

The above picture is confirmed when we simulate ac- 
tual multi-layer growth. The cut through a crystal de- 
picted in Figure [j] shows the crystal converging to the 
stable thickness within 5-10 layers and then continuing 
to grow at that thickness. 

We can understand better the reasons for this be- 
haviour by examining example polymer configurations 
and the probability distributions of the stem length for 
the growth of a single layer at a number of surface thick- 
nesses (Figure ||). l m i n places one constraint on the stem 
length; only a small fraction of the stems can be shorter 
than / m ; n if the layer is to be thermodynamically sta- 
ble. The boundary of the growth face places the sec- 
ond constraint on the stem length; it is energetically un- 
favourable for the polymer to extend beyond the edges 
of the underlying surface. Even in the absence of this 
constraint, i.e. on the infinite surface, the stem length 
remains finite because at every step there is always a fi- 



nite probability that a fold will be formed jl(|[ll| , and so 
the probability distribution decays to zero at large stem 
lengths (Figure pi) . 

The effect of the finite thickness of the growth face can 
be clearly seen in the growth of a layer on a 50-unit thick 
surface (Figure ||b); the probability of the stem length be- 
ing greater than the surface thickness is much less than 
it being smaller. Therefore, the new layer is thinner than 
the underlying surface. Only when the surface thickness 
approaches l m i n does the probability distribution become 
symmetrical about the surface thickness and the thick- 
ness of the new layer become equal to the surface (Figure 
. This is the reason why the thickness at which stable 
growth occurs is close to l m [ n . 

The above scenario does not hold for all temperatures. 
The last point of the curve in Figure ||a is for a surface 19 
units thick; a new layer cannot grow on a thinner surface. 
However, there is no a priori reason why the thickness 
curve must end after the curve has crossed y — x. In- 
deed, for T > 3.2 efc -1 there is no thickness for which 
successive layers have the same thickness. For example, 
after the growth of two layers on a 50-unit thick surface 
the outer layer is ~29 units thick (Figure ^>); the crystal 
then stops growing because the outer layer is too thin. 
For these smaller supercoolings, as suggested in the en- 
tropic barrier model, the entropic rounding of the crystal 
profile inhibits growth. 

To overcome this entropic barrier requires a cooper- 
ative mechanism whereby a new layer takes advantage 
of (and then locks in) dynamic fluctuations in the outer 
layer to larger thickness. The presence of such fluctu- 
ations is shown in the probability distributions of Fig- 
ure ||; however, growth stops in our model because we 
attempt to grow a new layer on an outer layer that is 
static. As overcoming the barrier would be most rapid 
when the magnitude of the fluctuations is the minimum 
necessary, we expect that this mechanism would lead to 
the crystal continuing to grow with a thickness that cor- 
responds to the smallest for which a new layer can grow, 
e.g. at T = 3.375 efc -1 this is a thickness of 30 (Figure 
||b). This mechanism again leads to a thickness for the 
polymer crystal which is close to l m i n - 

In summary, our results present a new picture of the 
mechanisms that cause the lamellar thickness to be con- 
strained to a value just above l m i n - The free energetic 
costs of the polymer extending beyond the edges of the 
previous crystalline layer and of a stem being shorter 
than l m i n provide upper and lower constraints on the 
length of stems in a new layer. Their combined effect 
is to cause the crystal thickness to converge to a value 
close to l m i n where growth with constant thickness then 
occurs. 

This convergence of the thickness has been observed in 
experiments in which the thickness of growing polymer 
crystals adjusts to a change in temperature Q and in 
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FIG. 4. Polymer configurations and probability distributions of the stem length for a new crystalline layer grown at 
T — 2.75 ek" 1 on a surface of thickness: (a) oo (b) 50 and (c) 21. Each configuration contains 150 stems and a dashed line 
marks the boundaries of the underlying surface; the y-axis origin is at the centre of the surface. The dashed vertical line in 
the probability distributions in (b) and (c) is at the thickness of the underlying surface. 



which lamellae form by epitaxial crystallization onto 
extended-chain crystalline fibres [fl5| . It would be very in- 
teresting if atomic force microscopy ( AFM) could be used 
to probe the profiles of the steps on the lamellae that re- 
sult from temperature changes; from these profiles maps 
similar to Figure ||a could be constructed. AFM could 
be also used to study the profile of the crystal close to 
the growth face to examine whether entropic rounding of 
the crystal edge occurs. 

Our description of polymer crystallization differs sig- 
nificantly from that given by the surface nucleation ap- 
proach. The assumption that a crystal of any thickness 
can continue to grow at that thickness does not hold. 
Furthermore, although overcoming the free energy bar- 
rier to the nucleation of the first few stems in a new layer 
can be one of the slowest processes in the growth of a new 
layer, it does not play a role in determining the thickness 
at which stable growth can occur. By contrast, some of 
the insights of the entropic barrier model form part of 
the more complete picture we present here. 
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